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SUMMARY 

Net radiant heat flow between two infinite, parallel, tungsten 
plates was computed by summing the monochromatic energy exchange; the 
results are graphically presented as a function of the temperatures of 
the two surfaces. In general these fluxes range from approximately 8 
to 25 percent greater than the results of gray-body computations based 
on the same emissivity data. 

The selection of spectral emissivity data and the computational 
procedure are discussed. The present analytical procedure is so 
arranged that, as spectral emissivity data for a material become avail- 
able, these data can be readily introduced into the NASA data-reduction 
equipment, which has been programmed to compute the net heat flux for 
the particular geometry and basic assumptions cited in the text. 

Nongray-body computational techniques for determining radiant heat 
flux appear practical provided the combination of select spectral emis- 
sivity data and the proper mechanized data-reduction equipment are 
brought to bear on the problem. 


INTRODUCTION 

The radiant heat exchange between two infinite, parallel plates 
of opaque material can be readily computed by assuming that both mate- 
rials are gray bodies; that is, that their emissivity is dependent on 
surface temperature and not on the wavelength of the radiation. The 
more laborious method of summing the monochromatic energy exchange over 
all the wavelengths of significance, reported in reference 1, indi- 
cates that the gray -body method can yield sizable errors. For example, 
the gray-bodv method of computing the radiant heat transfer between 
infinite parallel plates of pure aluminum and Inconel over a wide 
range of temperature and temperature difference produced errors ranging 
from 2 to 29 percent. 

In the present paper the nongray-body heat flow between two tung- 
sten plates was computed on a data-processing machine using a finite 
difference formula. For any two given surface temperatures the net heat 



flux was determined by summing the net monochromatic heat flux over 
about 100 increments of wavelength. 

Tungsten was chosen for two reasons. Its radiant behavior has 
been studied more than that of any other material. It is a strong con- 
tender for use s.s cathode and anode surfaces in thermionic energy con- 
version devices. Such devices operate at hign temperatures., where the 
parasitic loss due to radiant heat exchange can have a pronounced 
effect on cycle efficiency. 

bd 

The hemispherical spectral emissivity data used in the heat-flux 
equation were obtained from several sources in an endeavor to arrive at 
the most accurate set of data obtainable at tnis time. Hemispherical - 
total -emissivity values were computed from these spectral data and are 
compared with the generally accepted values of total emissivity. The 
computed total-emissivity data were used to ootain gray-body flux rates, 
and the results are compared with the nongray-body results. 

The calculation procedure used in this report is described in some 
detail. The sources, and modifications thereto, of the tungsten emis- 
sivity data are presented along with a table of the hemispherical spec- 
tral emissivity data that were put into the heat-transfer calculations. 

A table and a figure show the net heat exchange for a. broad range of 
surface temperatures. 


ANALYSIS 

This section presents the basic assumptions included in the anal- 
ysis, the derivation of equations, the sources and modifications to the 
emissivity data, and the computational procedure. 


Assumptions and Conditions 

The calculations for nongray-body heat exchange between two in- 
finite, parallel, flat plates of opaque material incorporated the fol- 
lowing five assumptions or conditions: 

(1) The two plates are identical. 

(2) Each surface is at a constant tempeiature, surface 1 possessing 
a higher temperature than surface 2. 

(3) Each surface is a diffuse radiator. 
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(4) At a given temperature and wavelength the hemispherical spec- 
tral emissivity e \ ^ equals the hemispherical spectral absorptivity, 


(5) A vacuum exists between the two plates. 


Equations 

The basic equation (eq. 31-157 of ref. 2) for determining the 
nongray-body net heat exchange in watts per square centimeter is 


a 


whe re 

= 37,404 X -5 ( e 14, 387/XT _ x ) w/ (micron) (cm 2 ) 




in the appendix.) 


is the interchange factor. 


(Symbols are defined 


The equation written in finite difference form (assuming energy at 
wavelengths less than 0.2 micron as negligible and using the gray-body 
method at wavelengths greater than 20 microns) becomes 
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where a = 5.oG99 x 10 - 1 - 2 w/(sq cm)( 0 K^), AX =* 0.05 micron, and £ go rp 

the emissivity at 20 microns, is taken a.s the mean value of the emis- 
sivitv for the bundle of energy at wavelengths greater than 20 microns. 

The calculation of equation (2) requires that the spectral emis- 
sivity be provided at wavelengths of 0.2 to 20 microns in increments of 
0.05 micron for temperatures T]_ and Tg. 

The total hemispherical emissivity at the particular temperature 



When expressed in finite difference form, this equation becomes 


-*oT 4 (3) 


where AX = 0.05 micron. 

The equation used to determine the gray-body net-heat exchange is 



where the emissivity of the colder surface was evaluated at the geo- 
metric mean temperature 

T* = J T x T 

f 1 2 

according to the empirical method of reference 3, and was obtained 

from a temperature -emissivity plot of the results by equation (3). 

These gray-body data (eq. (4)) are used herein for comparison purposes 
only. 



t?d 
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Emissivity Data 

Hemispherical spectral emissivity data for tungsten were obtained 
as follows: 

Data at wavelengths less than 2.6 microns . - The normal spectral 
emissivity data obtained by DeVos (re f . 4) were selected for several 
reasons. Measurements were made over sizable ranges of temperature 
and wavelength. The quality of these data appeared very good. Emis- 
sivity values obtained from reference 4 data curves for temperatures 
of 1600° , 2200°, and 2800° K are shown in table I. The data were 
"corrected” in the following manner: 

(a) Reduced 2-| percent to correct for light scattering (ref. 5) 

(b) Changed from normal to hemispherical by means of figure 1 
(obtained from ref. 6) 

(c) Linearly interpolated or extrapolated to temperatures of 1000°, 
2000° , and 3000° K 

(d) Linearly extrapolated from 2000° and 1000° K to 0° K and from 
2000° and 3000° K to 4000° K (for reasons to be explained in the 
section Calculation Procedure). 

Item (a) was suggested by Larrabee, whose carefully controlled 
experiments were limited to the visible spectrum. His results (ref. 5) 
show good agreement with the results of reference 4. 

Converting the normal data to hemispherical data required only a 
small correction (fig. l), since most of the emissivity values (table i) 
were reasonably large. Nevertheless, two types of checks were made on 
the validity of the conversion curve. Reference 5 presents spectral 
emissivity data at angles of incidence from 0° to 90° for a wavelength 
of 0.665 micron. An algebraic integration of these data over a hemi- 
sphere yielded the data point shown in figure 1. The other check was 

to compute e 'v - and e. m using the self-consistent electromag- 
K £ a. , I , n 

netic wave equations in reference 6. Admittedly, both equations, while 
reasonably good at the longer wavelengths, yield inaccurate spectral 
emissivity values at the shorter infrared wavelengths of interest. 
Nevertheless, since the equations for ip e \,T,n were self- 

consistent, it was felt that they would be indicative of the validity 
of the conversion curve as applied to the spectral data for tungsten. 
Calculations were made at 2000° and 3000° K (over a range of wavelength 
from 1 to 2.5 microns). Results, not shown, were in excellent agree- 
ment with the curve. 

Since e varies appreciably with temperature at most wave- 

X 9 I 

lengths of interest, it is fortunate that a nearly linear relation 
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exists between e^ T and temperature (table i). Hence, linear inter- 
polation and extrapolation of these data to temperatures of 1000°, 

2000°, and 3000° K appear well justified. 

Data at wavelengths greater than 4 micron s. - In this region 
experimental values of spectral emissivity were not found in the lit- 
erature for broad bands of temperature and wavelength. Equations of 
Schmidt and Eckert (ref. 6) were considered to be as reliable as any 
other at long wavelengths for temperatures between 1000° and 3000° K. m 

The procedure was as follows: 

(a) Equations 4-84 and 4-85 of reference 5 were used to obtain the 
data for temperatures of 1000°, 2000°, and 3003° K (table II). The 
equations arei 

For 0 < (r T /x)l0 -2 < 0.5 

e \,T = °* 476 |/"(i»tA)10 _2 - 0.148 Q^)x 10 -2 (5a) 

For 0.5 < (r T /X)l0"2< 2.5 

e X,T 3 °-442 (np/x ) 10 -2 - 0.0995^) x |_0" 2 (5b) 

Data for the electrical resistivity r»p were obtained from references 
7 and 8. 

(b) Then these data were linearly extrapolated to 0° and 4000° K 
in the manner described for the short -wavelength data. 

Data at wavelengths between 2.6 and 4 mic uns . - The T values 

were selected by plotting the previously described two sets of data 
(fig. 2) and joining them by a smoothly faired curve. 

Surface properties . - To describe the surface properties of the 
metal represented by the emissivity data used :.n this report would be 
to toy somewhat with fancy. Yet some sort of a statement is called for 
if the heat-flux data are to be of engineering value. Data at the long 
wavelengths were for a "smooth" surface. Perusal of figure 3 will show 
that, at a temperature above 2000° K, more than half the energy exchange 
will be at the shorter wavelengths where data of reference 4 apply. 

DeVos painstakingly defined the pertinent properties of the rolled 
strip. Briefly stated, the tungsten was of high purity and was prop- 
erly outgased and then annealed to ensure that no systematic changes in 
emissivity would occur after the annealing process. The final crystal- 
line texture showed the crystal grains to be of the order of 30 microns 
and the surface to be predominately of the (lOO) plane. 
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Calculation Procedure 

The requirement of a small wavelength interval AX in order to 
obtain a good precision of results makes automatic data-processing 
equipment ° a practical necessity. The program used herein was so 
arranged that the equations to be solved were placed on. one input and 
the tungsten emissivity data were placed on another. With the compu- 
tation program arranged in this manner, little effort is required to 
change the spectral data input, and therefore the program. is useful for 
computing the heat flux of other materials as spectral emissivity data, 
become more plentiful. (The present program does not require that both 
surfaces be of the same material. ) 

Size of wavelength interval. - The selection of a wavelength in- 
terval AX equal to 0.05 micron was based on the Planckian energy dis- 
tribution (see fig. 3). Even with this relatively small interval, as 
much as 3 percent of the black-body energy at 3000° K is contained in 
an interval. However, the interval selected appeared satisfactory, 
since there was no noticeable data, scatter observed during the plotting 
of the results of equation ( 2 ). With an interval severalfold larger 
than the one used, precision of the results would have been compromised. 

Cutting off the orderly summation at a. wavelength of 20 microns 
was justified on the basis Of figure 3. At 1000° K the black-body 
energy in the tail (see eq. (2)) is about 1.5 percent of' the total 
energy. Furthermore, the spectral emissivity decreases with increasing 
wavelength. 

Temperature input . - The temperatures selected for computation of 
q by equation (2; were 


0 

1 — 1 

t 2 , °k 

4000 

3800, 

3600, 3400, . . 

. 600 

3800 

3600, 

3400, 3200, . . 

. 600 

3600 

3400, 

3200, 3000, . . 

. 600 

800 

600 




Equation (4), used to obtain the total hemispherical emissivity, 
was evaluated at temperatures of 1000°, 2000°, 3000°, and 4000° K. 

Emissivity input. - Tungsten emissivity data were obtained directly 
from tables I and II. These data, as introduced into the computing 
machine, are presented in table XXI. A linear interpolation procedure 
was used to obtain e ^ p at the particular values of wavelength and 

temperature needed in equations ( 2 ) and (3). 




The emissivity data were extrapolated to temperatures above 3000° K 
and :>e low 1000° K (tables I and II) in order to obtain flux data over a 
wider range of temperatures than would have been possible otherwise. At 
temperatures between 3000° K and the melting joint, the linear extrap- 
olation of the 2000° to 3000° K data is reasonably well justified. At 
these high temperatures (short wavelength region, i.e., table I data) 
the spectral emissivity varies almost linearly with temperature. At the 
other temperature extreme, the error in spectral emissivity associated 
with the linear extrapolation of the 2000° to 1000° K data into the 
region below about 800° K could be appreciable. Hence, flux calcula- 
tions were not made for values of Tj lower than 800° K. 

The selection of wavelengths for the data input, given in table III, 
followed no hard and fast rule; however, greatest weight was given to 
the emissivity data at the shorter wavelengths, since the majority of 
the radiant energy is emitted there (figs. 2 snd 3). 


RESULTS AND DISCUSSION 

The nongray-body heat -transfer data obtained using equation (2) are 
presented in their entirety in table IV. Also, they are shown in graph- 
ical form in figure 4. There the lines of constant temperature T-^ 

are shown as dashed lines at all values of T^ - T ? that are outside 

the temperature framework of the table. Data at temperatures above the 
melting point are useful only insofar as they can facilitate interpo- 
lation of Tj data, near the melting point, which is shown in the figure. 


Under the heading "Emissivity Input" the possibility of a rather 
large error in spectral emissivity for low values of temperature is 
discussed. It is primarily for this reason that the curves of constant 
temperature Tp in figure 4 are shown as dashed lines in the vicinity 

of Tg equal to zero. The computed data (shovn as the solid lines of 

constant T^, when plotted as shown, lead one to believe the flux is 


almost constant for all values of Tg less than 600° K. Very likely 
this is the case. But it should be pointed out that, while the energy 


{£ 


is relatively small. 


emitted by the To surface 

’ ‘ l 2 

the amount of energy the T 2 surface absorbs Ls determined in consid 
erable measure by its spectral absorptivity. As a rough approximation 

/ QD 

^X,Tp e X,T]^ e X,T 5 ^X. For l ar £ e values of 
Tj, the values of UjTg that have the greatest influence on the net 
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heat flux are those having short wavelengths (fig. 2). Unfortunately, 
the rn values in the low- temperature and short -wave length region 

were subject to a large temperature extrapolation (table i). This ex- 
ample illustrates the caution that must be exercised when one uses heat- 
flux data based on a rather limited amount of spectral emissivity data. 

In figure 5 hemispherical total emissivities computed from the hem- 
ispherical spectral data of table III by means of equation (3) are com- 
pared with the widely used data of reference 9. The comparison of 
emissivities (fig. 5) is not intended to show either set of data to be 
in error. Rather, it does show that both sets of data are in fair 
agreement and that gray-body heat-flux calculations made using either 
set of data would give similar results. 

A comparison of the gray- and nongray-body methods of computing heat 
flux is given In figure 6. The emissivities used in the gray-body cal- 
culations were obtained from the curve shown faired through the four 
data points given in figure 5. There is very little data scatter ob- 
served in figure 6. This attests to the desirability of using small 
increments of wavelength AX in making radiant -he at -flux calculations 
by means of finite difference formulas. In all cases, the nongray-body 
flux is shown to be greater than the gray-body flux. In general, the 
discrepancy increased as Tq was reduced and/or as Tq - T£ was in- 
creased. The minimum deviation occurred near the melting point and was 
approximately 8 percent. The maximum was about 25 percent. These re- 
sults can be compared with the results obtained in reference 1 for In- 
conel, whose spectral emissivity was shown to decrease linearly with 
increasing wavelength. For most cases reported, the gray-body flux data 
ranged from 20 to 25 percent lower than the nongray-body flux. 


SPECTRAL DATA ON OTHER MATERIALS 

Reasonably accurate spectral emissivity data on high -temperature 
materials, although still greatly limited in quantity, are currently 
becoming available. It is most desirable that experimental data be 
obtained over the range of wavelengths necessary for heat-flux calcula- 
tions. As an example, the data points in figure 2 are from an inves- 
tigation (ref. 10) wherein normal spectral emissivity measurements were 
reported for a continuous band of wavelengths stretching from 0.45 to 
12 microns. This reference shows that the variations of spectral emis- 
sivity with wavelength for tantalum and molybdenum are similar to that 
of tungsten. Hence, these two metals may exhibit a. behavior not too 
much unlike that shown for tungsten in figure 6. 

Spectral data such as these, when obtained over a wide range of 
temperature and for surfaces of w r ell defined composition, will permit 
the computation of nongray-body flux data that are of greater accuracy 
than those obtainable by the gray-body method. 
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SUMMARY OF RESULTS 

A comparison of the nongray- and gray-bocy methods of computing 
net heat flux between two infinite, parallel jlates of tungsten shows: 

1- Gray-body calculated flux was from approximately 8 to 25 percent 
too low. Errors similar to these have been observed elsewhere when the 
variation of spectral emissivity with wavelength is large. 

2. Since many metals have a spectral emissivity - wavelength rela- 
tion similar to that of tungsten, appreciable errors in the net heat 
flux betw T een surfaces of these metals can be expected when grav-body 
calculation procedures are used. 

3. Experimental apparatus has been perfected to the point where an 
ample amount of spectral emissivity data on many materials can be ob- 
tained with sufficient accuracy to justify using nongray-body techniques 
of computing net radiant heat flux. 

4. The present analytical procedure is sc arranged that, as spec- 
tral emissivity data for a material become available, these data can be 
readily introduced into the NASA data- reduction center which has been 
programmed to compute the net heat flux for tire particular geometry and 
basic assumptions cited in the text. 


CONCLUSION 

The practicality of using nongray-body ccmputat ional techniques to 
determine radiant heat flux appears most encouraging provided the combi- 
nation of select spectral emissivity data and the proper mechanized 
data-reduction equipment are brought to bear cn the problem. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 15, 1961 


LLZl- 
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APPENDIX - SYMBOLS 

Planck’s distribution law equation, w/ (micron) (sq cm) 
net heat-transfer rate for gray-body conditions, w/(sq cm) 
net heat-transfer rate for nongray-body conditions, w/(sq cm) 
electric resistivity, (microhm) (cm) 
surface temperature, °K 

V T 1 x T 2, °K 

hemispherical total absorptivity 
hemispherical spectral absorptivity 
hemispherical total emissivity 
hemispherical spectral emissivity 
normal spectral emissivity 
w ave le ngt h } microns 

Stefan -Boltzmann constant ^ 5.6699*10 ^ w/ (sq cm) (°K^) 


hotter surface 


cooler surface 
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TABLE II. - TUNGSTEN EMISSIVITY DATA BASED 
ON ELECTROMAGNETIC THEORY 


Wave- 
length , 

* 

Temperature y T , 

°K 




r 



microns 

1000 

2000 

■ 3000 

■ 

0 

(a) 

4000 

(a) 


Hemispherical 

spectral 

e X,T 

emissivity, 

4 

0.112 

0.163 

0.1995 

0.061 

0.2355 

5 

.101 

. 147 

.1815 

.055 

.216 

6 

.092 

. 136 

.169 

.0475 

.2025 

8 

.081 

. 120 

.149 

.0425 

. 179 

10 

.0725 

.1085 

.135 

.0365 

.1615 

12 

.0675 

. 100 

. 123 

.035 

.146 

16 

.059 

.0875 

.110 

.030 

.132 

20 

.052 

.077 

.098 

.027 

.119 


®Ext rapo lat e d . 
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TABLE III. - HEMISPHERICAL SPECTRAL 


EMISSIVTTY DATA FOR TUNGSTEN 


Wave- 

Temperature, T, 

°K 

length , 





0 

2000 

4000 

inj-d uxis 





Hemispherical spectral 


emissivity, 




e \,T 


0.2 

0.4996 

0.4632 

0.4218 

.3 

.4996 

.4632 

.4218 

.4 

.4965 

.4628 

.4344 

.5 

.4856 

.4518 

.4221 

.6 

.4731 

.4398 

.4100 

.7 

.4639 

.4296 

.3916 

.8 

.4652 

.4153 

.3745 

.9 

.4485 

.4000 

.3608 

1.0 

.4205 

.3824 

.3545 

1.1 

.3934 

.3656 

.3470 

1.25 

.3404 

.3404 

.3404 

1.4 

.2995 

.3182 

.3361 

1.6 

.2477 

.2910 

.3322 

1.8 

.2052 

.2684 

.3266 

2.0 

.1710 

.2492 

.3223 

2.2 

.1439 

.2321 

.3174 

2.5 

.1169 

.2141 

.3087 

3.0 

.0882 

.1890 

.2850 

3.5 

.0710 

.1730 

.2500 

4.0 

.0610 

.1630 

.2355 

5.0 

.0550 

.1470 

.2160 

6.0 

.0475 

. 1360 

.2025 

8.0 

.0425 

.1200 

.1790 

10.0 

.0365 

. 1085 

.1615 

12.0 

.0350 

.1000 

.1460 

16.0 

.0300 

.0875 

.1320 

20.0 

.0270 

.0770 

.1190 




o 

o 

a ) 


tO O tO 
CO O) CM 

m co lo 

lo r~- a? 


^ H (O O 

cm p- co oo 

rM co co CD 


H ^ Od) O 
W tO LTJ K) ^ 

O rH to co CO 

N tO tO tO tO 

CO CO CM LO- K) LO 

to ^ lO CM CO rH 

to to ^ Oj o o 

to LO CO (C t- O 

cn cm 

W to ^ O CO ^ to 

H ^ ^ LO O) O O 



O 

« 

o 

o 

CM 


CM 

r -i 


EH ! 

o 

•s 

o 

Q) '■ 



CM 



■p 


o3 

o 

u 

o 

<v 

CO 

ps 

CM 

JE: 


0J 


Em 

O 


o 


c.D a> r- 
O tH CD 


Mi f— I 

rH CO tO 


O > H c; 
^ cd to 


to i— < CD CO to 

CM 'M 4 LQ CO r- 


CO O O co 

CD LO CD LO 


CM rH LO LO to 
LO t"- CO CD O 


O O CO CO CD CM ^ 
CO CD CO CO CM tO rH 


CO tO O 
LO LO tO 


O to CM fH 

r- O cm to 


CD CD CD 
CM CM CM 


00 CD CD 

^ to o 


to to t o 

H 4 •'tf 4 H 4 


co cm cm o 

r-i tO ^ LO- 


H tO O tO ^ 
LO CD O ^ CD 


CO o CO 

CM i — ! CO 


LO LO ^ 
CO co CD 


CD H tO ^ ^ ^ ^ 

O 1 — 1 <H rH ( — ( f— • i — t 

r~ L i — t i — L i — t i — I i — I i — I 


CM O- CO H 4 CD 


tO CM ' jl rH r-- 


CD r- H CM CD 

<0 O CM tO tO 


LO CO O rH H rH O 
H 4 ^ LO LO LO LQ LO 


CO CD 
G3 O CO 


CM p- CD LO «H 


o- t" co cn 


to O CD 
H 4 CO O 


tO H O CO M 1 
tO LO CO r- CO 


CD CM H 4 LO L0 LO L0 
CC CD CD CD CD CD CD 


IQ 'M 4 

M 4 tO tO W 


CD LO CO H 4 (H 


CD CD CD CM O 


CM CO CM CM 
LO CD tO CO 


o 

O L0' ^ co ^ 


lo H 4 to a> r - 

CD O H N to 

rH CM CM CM CM 


CM CD H 4 CM O 


CM LO P- CC- CD CD CO 

H 4 ^ H 4 tJ 4 ^ -"si 4 nJ 4 

CM CM CM CM CM CM CM 


CM -M 4 CO H 4 H 4 

CO rH L0 CD CM 

r H i — f i — I CM 

o o o o o 

o o o o o 

CO CD ^ CM o 

to tO to to to 


CL 1 CO i — L CM O 

^ to co CD o 

CM CM CM CM to 

o o o o o 

o o o o o 

CD (LO ^ 4 CM o 

CO CM CM CM CM 


LO CD H CM CM CM CM 

O O rH rH r-t rH rH 

tO to to to to to to 

o o o o o o o 

o o o o o o o 

CO CO Tf 4 cm o CO to 


-1277 




Ratio of hemispherical to normal emissivity 


r- 

r- 



0 ,1 .2 

Normal emissivity 


Figure 1. - Ratio of hemispherical to normal emissivity of 
metallic surfaces. Suitable for either total or spectral 
emissivities of tungsten. 





Figure 2. - Spectral emisaivity data for tungsten. 
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Figure 3. - Spectral distribution of black-body energy 
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Temperature, K 


Figure 5. - Comparison of hemi spherical total emis 
sivity data for tungsten. 




Ratio of gray-body to nongray-body heat 
exchange, q e /q ne 
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0 BOO 1600 2400 3200 

Temperature difference between hotter and cooler 


surface, - T?, °K 

Figure 6. - Comparison of gray- and nongray-body 
methods for computing net heat flux between two 
infinite, parallel plates of tungsten. 
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